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Abstract 

We derive the semiclassical limit of the coherent state propagator for systems with two degrees 
of freedom of which one degree of freedom is canonical and the other a spin. Systems in this 
category include those involving spin-orbit interactions and the Jaynes-Cummings model in which 
a single electromagnetic mode interacts with many independent two-level atoms. We construct a 
path integral representation for the propagator of such systems and derive its semiclassical limit. 
As special cases we consider separable systems, the limit of very large spins and the case of spin 
1/2. 

PACS numbers: 
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I. INTRODUCTION 



The spin-orbit interaction plays important roles in many areas of physics, from atomic 
physics to condensed matter. The quantum description of systems with such interactions 
requires the use of Hilbert spaces which are the direct product of the orbital space (for which 
the coordinate or the momentum eigenstates form a basis) times the intrinsic space of the 
spin. 

In quantum optics a similar situation arises in the study of the interaction between atoms 
and electromagnetic modes in a cavity. When only two states of the atoms are relevant, as 
for example in the ammonia maser, these two states can be formally represented by a spin 
1/2. The state of a set of N such atoms can be hkewise represented by the states of an 
angular momentum of larger magnitude. The Hamiltonian describing their interaction with 
a single electromagnetic mode of a cavity can therefore be written in terms of the operators 
a and a^, which annihilate and create excitations of the quantized electromagnetic mode, 
and Jz, J+ and J_, of the angular momentum. 

The semiclassical behavior of such systems has drawn attention for quite a long time. One 
natural representation for the study of this hmit is that of coherent states. The semiclassical 
hmit of the coherent state propagator for both the Weyl and the SU(2) group has already 
been studied in detail. The purpose of this paper is to derive the semiclassical limit of the 
coherent state propagator for general systems with two degrees of freedom in which one 
degree is canonical and the other a spin. 

The quantum propagator K{b"*,b',T) = (&"|e~'^^/^|6') represents the probability am- 
phtude that the initial state \b') be measured as after a time T, when evolved by the 
Hamiltonian H. The propagator is the essential ingredient for quantum dynamical calcula- 
tions and it is also fundamental in the study of the quantum-classical connection. Semiclas- 
sical approximations of the propagator in the coordinates and momentum representations 
were studied initially by Feynman [1] and later by many others [2]. Semiclassical formulae 
for the propagator in the basis of coherent states appeared for the first time in the works 
of Klauder [3, 4]. Although these papers have treated the propagator for both canonical 
coherent states \z) and spin coherent states not much work has been done on these two 
basis simultaneously. The development of the semiclassical theory occurred independently 
on each representation. 



2 



For canonical coherent states, Weissman [5, 6] re-derived the results of Klauder using the 
semiclassical theory of Miller [7] . In Weissman's work, and also in the original Klauder's pa- 
pers, the fluctuations around the critical trajectory have not been accurately performed, and 
a 'phase' factor was missed. In spite of this, a first numerical application of the semiclassical 
formula was performed by Adachi [8] for an two-dimensional chaotic map, obtaining good 
agreement with exact quantum results. The correct evaluation of the second order fluctu- 
ations appeared in the works of Baranger and Aguiar [9], Xavier and Aguiar [10, 11, 12] 
and, independently, by Kochetov [13]. However, a detailed derivation of the semiclassical 
coherent state propagator for one-dimensional systems has appeared only later in [14]. Ex- 
tensions of the formula for two dimensional degrees of freedom and applications to chaotic 
systems have been performed in Ref. [15]. 

There are two differences between the semiclassical formula of Baranger et al and that 
of Klauder and Weissman. First, there is the extra 'phase' in the new formula, which is 
usually complex and, in fact, is a signature of the basis of coherent states. It is related to 
the overcompletness of the basis, since changing the resolution of the unity also changes 
this term [14, 16]. A second difference consists in replacing the Weyl symbol H of the 
Hamiltonian operator H by the average Hz = {AH\z). This implies that the classical 
trajectories entering in the formula are subjected to instead of H. The dynamics with 
H actually appears naturally in the work by Klauder, but it was changed back to H to 
be consistent with the lack of the extra exponential factor of the formula. As discussed in 
Ref. [14], both these changes are essential to get good agreement between quantum and 
semiclassical results. We note that the semiclassical coherent state propagator also presents 
singularities and discontinuities due to phase space caustics and the Stokes phenomenon 
[8, 15, 17, 18, 19]. 

Approximations for spin coherent states have also been studied by Kuratsuji and Suzuki 
[20] . The basic difference between their approach and Klauder's is that the later represents 
the classical trajectories in a Bloch sphere by means of angle variables while Kuratsuji and 
Suzuki represents the same dynamics in terms of an stereographic projection of the Bloch 
sphere on a complex plane. As is the case of canonical coherent states, a detailed derivation 
in the spin coherent state propagator only appeared later with Solari [21], where features 
similar to those appearing in the canonical case were found: an extra exponential term, and 
the underlying classical dynamics governed by the average Hamiltonian Hs = {s\H\s). Vieira 
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and Sacramento [22] and Kochetov [23] have also derived the same formula independently. 
Yet another detailed derivation has also been presented by Stone et al [24] , focusing on the 
importance of the extra term (see also Ref. [25]), which has received the name of Solari- 
Kochetov phase. The singularities and discontinuities of the semiclassical spin propagator 
are discussed in Refs. [26, 27]. 

In spite of the extensive work on the semiclassical theory for the canonical and spin 
representations, there are very few results considering the two bases simultaneously. Alscher 
and Grabert [28] have calculated the semiclassical coherent state propagator for the spin 
1/2 Jaynes-Cummings model showing that it corresponds to the exact quantum result. 
Pletyukhov ct al [29, 30] derived a semiclassical trace formula for systems with spin, but did 
not consider the semiclassical propagator itself in detail. An overview of the semiclassical 
approaches for spin-orbit interactions and associated trace formulae was recently published 
by Amann and Brack [31], and a general semiclassical theory for Hamiltonians which are 
linear in spin operators has also been formulated [29, 30, 32, 33]. 

In this paper we derive a semiclassical formula for general canonical-spin systems using 
path integrals. We show that our formula agrees with the previous results [14, 21, 22, 23] 
in the separable case (no spin-orbit interaction) and that it reduces to the 2-D canonical 
coherent state propagator in the limit of very large spins. Finally we discuss the limit of 
validity of the approximation for spin-| systems. 

II. BASIC DEFINITIONS 

The propagator in the canonical-spin representation is given by 

K {z"*, s"*, z', s', T) = {z'\ s"\e-'"'^"'\z', s'), (1) 

where \z,s) = \z) (g) \s) is the product of a canonical coherent state \z) and a spin coherent 
state While \z) is defined in the 'particle' Hilbert space, |s) is defined in the {2j + 1)- 
dimensional space of an angular momentum j. These states can be written as 

\z) = e^^^-il^l^|0) and |.) = -^^\ -j), (2) 

(1 -I- isfY 

where z and s are complex numbers, is the canonical creator operator, j+ is the raising 
spin operator, |0) is the harmonic oscillator vacuum state and | —j) is the extremal eigenstate 
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{zi\z2) =exp<^ -T^Pil +ZiZ2 - -\Z2\ > and (si|s2) = 7- , (3) 



of J3 with eingenvalue —j. The coherent states are non-orthogonal with 

{i + \si\^y ii + \s2\'y 

and satisfy 

y ^|^)(^| = 1(.) and J -^—^\s){s\ ^ 1(,), (4) 

where x and y are the real and imaginary parts of z and X and F the real and imaginary 
parts of s. Throughout this paper we shall use lower case letters to refer to the canonical 
variables and corresponding upper case letters to refer to the spin variables. 

Finally, the complex number z labelling the canonical coherent state can be written 
explicitly in terms of position and momentum variables as 



V2 \b c. 

where q and p are the average values of the position and momentum operators, respectively, 
and the variances b and c satisfy be — h. 



III. PATH INTEGRAL FORMULATION 

In this section we construct a path integral representation for the quantum propagator (1). 
As usual we divide the time T into small intervals of size e = T/N and insert resolutions 
of unity between each propagation step. We obtain 

A-(."-./'-..'./.r) ^ .i./n{(^)^«p} 

7V-1 

l[{zk+i,Sk+i\e-'"'/''\zk,Sk), (6) 



7V-1 

X 

k=0 



where we define l^o, Sq) = \z', s') and {zn, sn\ = {z", s"\. In the limit e — > the infinitesimal 
propagators can be written as 

{zk+i, Sk+i\e~'^^'^^^\zk, Sk) = {zk+i, Sk+i\zk, Sk) exp i^--eHk+i/2^ , (7) 

where 

fv _ {zk+i, Sk+i\H\zk, Sk) , s 

J^k+1/2 = —, \ r— 1,0; 
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and 



(1 + 



2j 



1, ,2 * 1, ,2 

(1 + (1 + IsfcPy I 2 ^ 



With these considerations the propagator can be written as 

N-l nN~l 



K{z"*,s"*,z',s',T) = \im 



where 



with 



2j + l 



Yl {dxkdykdXkdYk} , 



k=l 



k=0 



(9) 



(10) 



(11) 



and 



■ N-l 
fe=0 



ihj X In 



fe=0 



(i + 4+i«fe)^ 



2i;i In 



(i + 4«ik) 



(12) 



(13) 



Eq. (10) is a discretized path integral representation of the propagator (1). In the following 
sections we shall consider the formal semi classical limit ^ — > and j ^ 00, keeping the 
product jh = J constant. 



IV. THE SEMICLASSICAL LIMIT 



In the semiclassical limit, the integrals in Eq. (10) can be performed by the saddle point 
method. The method consists basically in approximating the exponent F by a quadratic 
form and performing the resulting Gaussian integrals. The quadratic form is obtained by 
expanding F around its critical points. In the next subsections we shall: (a) find the critical 
points of F, and therefore the critical path; (b) compute F at the critical path; (c) expand 
F to second order around the critical path and compute the Gaussian integrals. 



A. The Critical Path 



We begin by looking for critical points of F. They satisfy the equations 
OF OF OF OF 



dzjYi 9zjj^ dsjYi 



0, m = 1,...,A^- 1. (14) 
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As Hk+i/2 — -f^fe+i/2(-2fe+i, sl^i, Zk, Sk), these equations can be explicitly written as 

ie dHjn+i/2 



ie dH, 



m-1/2 



n dz^ 



^m— 1 ) 



and 



ie dHm-l/2 . I Sm Sm-l 



-J 



(15) 



(16) 



It is important to emphasize that, because m = 1, . . . , iV — 1, the variables Zq, Sq, z^, sn 
do not enter in Eqs. (15) and (16): the critical path, defined by the set of critical points, 
depends only on zo, sq, z*^ and s^, and not on Zq, Sq, z^-, sm- 
In the limit e — > Eqs. (15) become 



i 9H _ , i dH . 

Ilk-' ""^ la? = -' 

where H = {z, s\H\z, s) . In terms of q and p (see Eq.(5)) this corresponds to the usual 
Hamilton's equations 

dH . ^ dH . 

— = , and — = (18) 

For Eqs. (16) the calculation is slightly more involved but the result is also very simple. 
We obtain, in the limit e — > and — > with j — J/h, 

dH 2iJs* , dH 2iJs , , 

— — = o and — — = o- 19 

ds (l + s*s)' ds* (l + s*s)' ^ ' 

The trajectories described by Eqs. (17) and (19) define the critical path of the Feynman 
integral (10). Nevertheless, these trajectories must satisfy the boundary conditions ^(0) = z\ 
s(0) = s', ^*(r) = z"* and s*(r) = s"\ as can be seen from Eqs. (15) and (16). As 
discussed in detail in [14], these trajectories are usually complex and the variables z and z* 
are not generally the complex conjugate of each other, the same happening between s and 
s*. Therefore, it is convenient to set a new notation 

z^u, s^U, z* and s* = V. (20) 
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In terms of these new variables, the equations of motion (17) and (19) become 
idH _ . idH __. idH _ 2jV i dH _ -2jU 

n~d^~'"' h~d^~~'^' ndu ~ + nW ~ {i + uvf ^ ' 

with boundary conditions 

m(0) = z', t/(0) = s', v{T) = z"\ ViT) = s"\ (22) 

Since z'* , s'*, z" and s" do not appear in the boundary conditions, the value of u{T), U{T), 
v{0) and V^(0) arc determined by the integration of Eqs. (21). From now on we shall use 
this new notation to refer to the complex classical trajectories. The discrete variables z^, 
Sjn, and shall also be replaced by u'^, U^, v'^ and V^, respectively. 



B. The Complex Action 



The function F appearing in Eq. (11) can now be calculated at the classical trajectory. 
We use a bar over the variables to indicate that they are calculated at the critical path. We 
have 

. N-l 



k=0 
. N-l 



(1 + V'U'') (1 + 



N-l 



k=0 



k=l 



II*, II II* -0-0 -N-N\ -1 

— ' z z — u V — u V ) — J In 



{I + s' s'*){l + s" s"*) 



[1 + WV^){1 + U^V^)_ 

As usual we have replaced z'* , s'*, z", s" by v^,V^, , in the first line and corrected for 
this in the last line. Taking the hmit e — > we find, after some algebra, 

N-l 



F = ^S{z"\ s"\ z\ s', T) - A - 2 ^ In (1 + V^U^) , 



(23) 



where S{z"*, s"*, z',s', T) is the complex action 



■T 



S{z"*,s"*,z',s',T)^ / { — {ilv - du) - ihj 



UV -VU 
1 + UV 



H\dt 



^ {u'v' + u"v") - ihj In [(1 + U'V'){1 + U"V")] 



(24) 
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and 

A = ^ (l^f + WT) + J In [(1 + \sf){l + Wf)] (25) 

is a 'normalization term'. The limit e — has not been taken on last term of Eq.(23) 
because this term is going to get cancelled later on. 
It can be checked that 

dS -2ihjV' dS -2ihjU" 



dU' l + U'V dV" 1 + V"U"' 

and 

dS 



(27) 



where a single (double) prime to denotes initial (final) time. 



C. Fluctuations around the Critical Path 

In the scmiclassical limit, the only relevant points in the path integral of Eq. (10) are the 
saddle points. In the present case, they define trajectories governed by Eqs. (21) [or by their 
discretized forms, Eqs. (15) and (16)]. The exponent F in Eq. (10) has to be integrated 
over the intermediate points x = {u^, U^, v^, V^, . . . , u^~^, U^~^,v^~^, yN-i^r^ Expanding 
F around the critical trajectory x we obtain 

F(x) = F(x) - -5x^[-52F(x)]5x, (29) 

where 5x = x — 5c. The matrix [— (5^F(x)] contains the second derivatives of F calculated 
with the critical trajectory. Substituting this result in Eq. (10) and considering the jacobian 

{dxkdvkdXkdYk} ^ ~ {dzkdzldskdsl} = {du'dv'dU'dV'} , 

we find 

K = lim I (^3^) J n {d[Su^]d[5v'']d[5U'']d[dV^]} e-i'5-^[-'5'-P(*)l'5x| ^ (39) 
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where -F(x) is given by Eq. (23). The matrix [— 5^F(x)] is written as 

/ n-iN-i 111 
rill 



iN-l 

'-21 

1^21 ''-22 

1 





1 


























-2 
















njN-2 
^33 


njN- 


-2 


^13 - 1 


njN- 

rT-23 


-2 








njN-2 


niN- 
''•44 


-2 


njN-2 
/t4i 


njN~2 
/T24 


j^N-2 








r ^ - 1 


^41 


-2 


-T/iV-2 
/til 


^21 


-2 


1 





a/Af-2 
rt23 


'i/JV-2 

ri42 


^^-2 


njN-2 
^21 


niN- 
1^22 


-2 





^JV-3 



v ■■ 

where 



2(J + 1) 
(1 + V'^U'^Y 



2j 



(1 + 



(31) 



and = d T-Lm+i/2/ dxT^ ^ ■ In this last definition we are using the compact notation 



Xi = X2 = U, X3 = V and X4 = V^, 



(32) 



In addition, rrii {rrij) equals to m + 1 when i (j) is equal to 3 or 4, and equals to m when 
i (j) is equal to 1 or 2. In order to cancel the last term in Eq.(23) and also the factor 
(2j + in the pre-factor of Eq.(30) we change the variables associated with the spin by 

the transformation SU"" = dU"'/VB^ and SV"" = SV'^/y/B^. Eq. (30) then becomes 



i^ = e^'5(x)-Aiiin 



I V 47r2 



N-l 



k=l 



<5x^[-<52f(x)]<5x 



,(33) 



where iS(x) and A are given by Eqs. (24) and (25) respectively and 5x 
{Su\ SU\ Sv\ SV\ Su^-\ SU^-\ Sv^-\ SV^-^ and [-S^F{5t)] is the matrix 



- 5^F{5t) 







v 



(34) 
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where 



11 



21 



and 



1 







nl^^ - 1 





1 



1 








1 


"33 




















lo/fc-l 
'T-23 






o\ 




0/ 



The matrix (R"^)^ is the transpose of R"* and = l/y/B^- 

The Gaussian integral in Eq. (33) can be solved immediately and the result is 



N-l . 

Yl |dMdMd[(5i/'=]d[5\/^]} exp I --Si^[-S^F{ft)]S: 

k=l ^ 



(271)4(^-1) 
det[-(52F(x) 



(35) 



The evaluation of the fluctuation determinant is the most lengthy step of the whole calcu- 
lation. Here we shall only briefly describe the main steps of the calculation. 

We call the determinant of [— 5^F(x)] and we solve it by the Laplace method of 
reducing it to smaller determinants. In this process we are lead to define 5 auxiliary ma- 
trices whose determinants, together with A^, form a closed set of six discrete recurrence 
relations. The five determinants are called A^ and A^, for i, j = 1 or 2. The former is the 
determinant of [— 5^F(x)] without the first two lines and columns, while A^ is the determi- 
nant of [— 5^F(x)] without the first, second and i-th hues and without the first, second and 
j-th columns. Taking the limit of ^ 00 in the mentioned set of relations we obtain the 
following set of linear differential equations: 





( ° 


— H22 


-Hn 


— H21 


— H21 





\ 




H44 


— 2H24 





— H41 


Hn 


Hn 




i 


H33 





-2Hi3 


H23 


— H23 


— H22 




h 


H43 


— H23 


— H41 


-H+ 





H21 






H43 


— H23 


— H41 





-H+ 


H21 








H33 


H44 


— H43 


— H43 


-2H+ 


/ 
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D 



(36) 



where D^(i) = Aii(i), A.22{t), Ai2(t), ^2i{t), ^o{t)) is a vector containing the original 

determinant A{t) and the five auxihary ones. Hij are elements of the matrix 



H = 



dudu 



d 



dudU 



dudv 



d 



d 



dUdu 

dvdu 



d 



dVdu 



d^ 



d^ 



dUdU 



+ 



2V dH 



l+UV dU 



d 



dvdU 



dUdV 



l+UV 



d 



d 



dUdv 

dvdv 
dVdv 



d2 



d^ 



dUdV 



dudv 
V 



+ 



dH , rrdH 

av~'~^ au 
l+UV 



d 



dvdV 



+ 



2U dH 



dVdV ' l+UV dV 



) 



calculated with the classical trajectory where (/=(! + UV)I \f22 and H-t = H13 ± H24. 
Setting D' = De~^/^+<^* we obtain the more symmetric form 



n 



(37) 



with 



A 



H44 
H33 
H43 



-H22 
H_ 



-H23 



Hn 


-H 
-H41 



H43 — H23 — H41 
\ H33 H44 



-H21 
-H41 

-H23 



-H43 



-H21 

-H41 

-H23 



-H43 



\ 

Hn 
— H22 
H21 
H21 



(38) 



The equations for D' are intimately related to the linearized equations of motion around 
the critical trajectory. To see this we go back to our notation as is Eqs. (32). Setting the 
small displacements 5xi(t) = Xi(^) — Xi(^) around the critical trajectory and defining 



ei(t) = bu{t), 6(t) 



[1 + U{t)V{t)) 



/2j 



SU{t), 



(39) 



(40) 



[l + Uit)Vit)) 

we can construct an anti-symmetric tensor Tik{t) = C,i{t)C,l.{t) — ^i{t)^kit), where C,i{t) 
and 4'(^) ^-re independent displacements from the critical trajectory. The tensor T has 
six independent components, which can be arranged into a new vector defined by — 
{Tu{t),T23{t),T4i{t),Ti3{t),T42{t),Ti2{t)) whose equation of motion is exactly T = iAT /%. 
Putting things together we find that 



det[-(52F(x)] = A(T) = r34(T) exp 



T 



(41) 
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Since T34 is related to the linearized motion around the critical trajectory, it can be 
easily written in terms of the tangent matrix or in terms of second derivatives of the action. 
Working out the details we find 



A{T) 



(42) 



(1 + U{T)V{T)) 

where M},h is the lower right 2 by 2 block of the tangent matrix in the coordinates xf- 

( Mn(T) Mi2(T) Mi3(T) M^,{T) \ 
M2l(r) M22(T) M23(r) M24(r) 

M3l(T) M32(T) M33(r) M34(T) 

Vm4i(t) M42(t) M43(r) M44(r)y 

Differentiating both sides of Eqs. (26) and (27) and conveniently re-arranging the terms we 
identify 



M(T) 



(43) 



M, 



ih 



bb 



du'dv" dU'dV" du'dV" dU'dv" \ gu'dv" \ 1+U'V' ) 



dU'dV" '^J \ \+U'V') du'dV" 
1 2 'd^S 



du'dv" 



(44) 



D. The final formula 



Collecting all the results for the exponent and pre-factor, the final formula for the semi- 
classical limit of the canonical-spin coherent state propagator becomes 

.1/2 



K{z"\s"\z\s\T) 

where 

S{z"\s"\z\s',T) = 



1 + U"V'' 



1 + U'V / det M 



1 



exp<{ -(5 + ^) - A 



^ ^ in ,. . , . I uv -vu , , , 

uv - vu) - ifij I ^ _^ I - Hy dt 





- {u'v' + u"v") - ifij In [(1 + U'V'){1 + U"V")] , 
2 



g{z"*,s"*,z',s',T) 



1 

2 JO 

1 
2 



dvdu 



H.dt. 



d {l + VUfdH d {1 + VUfdH 



dV 2j dU dU 2j 



dV 



(45) 



(46) 

dt 
(47) 



13 



and 

A = ^ (l^f + Wf) + i In [(1 + |.f )(1 + |.'f )] . (48) 

Alternatively, the pre-factor can be written explicitly in terms of derivatives of the action 
acording to Eq. (44), 

-/ l + U"V" \ 1 
\l + U'V' J det Mbb 

where 

au'av" du'dv" 

a^s a^s 
au'dv" au'av" 

All quantities are to be calculated at the stationary trajectory (we have removed the bar on 
top these quantities to simplify the notation). 



[1 + U"V"){1 + U'V') 
2i 



[det E] 



(49) 




V. SIMPLE APPLICATIONS 



In this section we shall apply the semiclassical formula Eq.(45) to three simple situations: 
(a) non-interacting spin and field operators; (b) the limit of very large spins and; (c) the case 
of spin 1/2. In each case we shall see that the propagator obtained from Eq.(45) reduces to 
well known results. 

A. Non-interacting Hamiltonian 

If the Hamiltonian can be separated into H = Hz + Hg, where Hz = Hz{a,a^) and 
Hg = Hs{J+, J-, Jz), then H = {z,s\H\z,s) = Hz + Hg, where Hz = {z\Hz\z) = Hz{z*,z) 
and Hs = = Hs{s*, s). Therefore, the complex action of Eq. (24) takes the form 

S{z"\ s"% z\ s', T) = Sz{z"\ z\ T) + Ss{s"\ s\ T), (51) 

where 

S,{z"*, z', T) = /(f {f {iiv - vu) - Hz] dt - {u'v' + u"v") , 

(52) 

Ss{s"*, s', T) = [-ihj (^^) -Hs}dt- ihj In [(1 + U'V'){1 + U"V")] . 
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Moreover, the term Q of Eq. (47) becomes 



(53) 



where 



gs{s"\s\T) = \j^ 

Finally, the det E simphfies to 

detE = — 



+ 



dV 2j dU ' dU 2j aV 



dt. 



du'dv" dU'dV" 



Therefore, for non-interacting systems, Eq. (45) amounts to 



K {z"\ s"\ z\ s\ T) = {z"\ z\ T) x K, {s"\ s\ T) , 



(54) 



(55) 



(56) 



where 



K,{z"\z\T) 



h du'dv" 



i {1 + U"V"){1 + U'V') d^S, 



n 



2j 



MSs+gs)-jln[(l+\s'n{l+\s"n] (rrr^ 



dU'dV- 



are exactly the semiclassical propagators known in the literature for the Weyl and SU(2) 
groups respectvely (see, for example, Refs. [14] and [24]). 



B. The limit of large spin 

Following Perelomov [34], wc set s = w/y/2j, J+ = \/2j^^ l^t j — >■ oo. In this limit 
the spin coherent states transform into canonical coherent states: 

In this case, discarding terms smaller than we obtain 
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1 + ss* 2^ 

J In [(1 + s's'*){l + s"s"*)] ^ {w'w'* + w'V'*) , (59) 
(9(1 + ss*)2 9^ ^ 9(1 + ss*f dH 2 



5s* 2j ds ds 2j ds* dwdw* 

and 



(l+s" s"*){l+s' s'*) j n dz'dz"* h dz'ds"* _^ [ h dz'dz"* h dz'ds"* I /gg\ 
\ nds'dz"* nds'ds"* / \ hdw'dz"* h dw'dw"* 

Equations (59) to (60) applied to Eqs. (45) and (49) produces the two-dimensional semi- 
classical propagator for canonical coherent states [15]. 



C. Spin-i Systems 

The semiclassical approximation developed in sections II to IV employed explicitly the 
limit j ^ oo. In this subsection we discuss the application of Eq. (45) to spin-| systems, 
whose general Hamiltonian is given by 

H ^Ho + Hs = Ho{a, a^) + M ■ C(a, a^). (61) 

In this case the classical Hamiltonian reads 

H{u,v,U,V) = {z\Ho\z) +h{s\s\s) ■ {z\C\z) = Ho{u,v) + -H,{u,v,U,V), (62) 



2 



where 



U + V V -U 1- uv 

H,{u, V, U, V) = C,{u, v) -^-^ - iC,{u, v) ^—^ - Cs{u, v)^—^. (63) 

and {z\C\z) = {Ci{u,v),C2{u,v),Cs{u,v)). 

The equations of motion are given explicitly by 

* = t£(^o + |^.), 



V = I [(Ci + iC^) - (Ci - iC2)V' + 2CsV] , 
if = I [(Ci + iC2)U^ - (Ci - ^C2) - 2CsU] . 



(64) 
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In the limit of small h we can drop the terms ^Hs on the first two equations and decouple 
u and V from the spin variables U and V. These, on the other hand, describe the precession of 
the spin in the external field C{u, v) generated by the orbital motion. In this approximation 

the orbital part of the action also separates from the total action and the semiclassical 
propagator can be written as 

K {z"% s"*, z', s', T) = {z"\ z\ T) x T) , (65) 

where is the one-dimensional canonical propagator and i^^s can be written as [23] 

^s[u,v]i^ ,s,T) = (i + |g//|2)(i^|g/|2) ■ (66) 

The coefficients a{t) and b{t) are obtained from the differential equation 

dW i 

^ = -1<7 • C{t)W{t) (67) 

where 

^ ^ f a{t) h{t) \ ^ ^ 

Wit) = ' ' ' M , (68) 
\-h*{t) a*it) J 

a are the Pauli matrices and W{0) = 1. Since Eq.(66) is the exact propagator for a spin in 
an external field, Eq.(65) can also be derived directly from the path integral approach 

K {z"\ s"\ z\ s', T) = J (."*, T) et^-["'^'^l (69) 

where the steepest descent approximation is performed only in the orbital action F^q (which 
contains only Hq). The spin propagator Kg is viewed as a slow varying pre-factor and is 
simply calculated at the stationary trajectory [32, 33]. This shows that the semiclassical 
formula Eq. (45) can also be used for systems with spin j = 1/2, in spite of the large spin 
limit considered in its derivation. 
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